Electrostatic Contribution to Twist Rigidity of DNA 
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The electrostatic contribution to twist rigidity of DNA is studied, and it is shown that the 
Coulomb self-energy of the double-helical sugar-phosphate backbone contributes considerably to 
twist rigidity of DNA — the electrostatic twist rigidity of DNA is found as C" c i cc ~ 5 nm, which makes 
up about 7% of its total twist rigidity (Cdna ~ 75 nm). The electrostatic twist rigidity is found, 
however, to only weakly depend on the salt concentration, because of a competition between two 
different screening mechanisms: (1) Debye screening by the salt ions in the bulk, and (2) structural 
screening by the periodic charge distribution along the backbone of the helical polyelectrolyte. It is 
found that depending on the parameters, the electrostatic contribution could stabilize or destabilize 
the structure of a helical polyelectrolyte. 

Pacs numbers: 87.15.-v, 36.20.-r, 61.41. +e 
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I. INTRODUCTION 

Genetic information in living cells is carried in the 
double-helical linear sequence of nucleotides in DNA. The 
DNA double-helix can be found in several forms that 
differ from each other in the geometrical characteristics 
such as diameter and handedness. Under normal phys- 
iological conditions, DNA adopts the B-form, in which 
it consists of two helically twisted sugar-phosphate back- 
bones with a diameter 2.4 nm, which are stuffed with 
base pairs and are located asymmetrically with respect 
to each other as characterized by the presence of major 
and minor grooves. The helix is right-handed with 10 
base pairs per turn, and the pitch of the helix is 3.4 nm 

[!]• ^ 

It is well known that above pH 1 each phosphate group 
in DNA has a negative charge [2], which renders the 
polymer stiff due to the electrostatic repulsion between 
these groups. The presence of neutralizing counterions 
and salt in the solvent screens the electrostatic repulsion, 
thereby leading to an effective way of controlling the stiff- 
ness of polyelectrolytes via the ionic strength of the so- 
lution. To account for the electrostatic stiffening, Odijk 
[3], and Skolnick and Fixman [4] have adopted an effec- 
tive wormlike chain description for the bending elasticity 
of stiff polyelectrolytes, and have calculated the correc- 
tion to the persistence length due to the electrostatic in- 
teractions. The so-called electrostatic persistence length 
is found to be proportional to the square of the Debye 
screening length, implying that the stiffness of polyelec- 
trolytes such as DNA should be very sensitive to salt con- 
centration [5] . While there are experiments that measure 
the electrostatic contribution to the bending rigidity of 
DNA in various salt concentrations [6,7], it is generally 
believed that changing the ionic strength has no signif- 
icant effect on the rigidity of DNA under most physio- 
logically relevant conditions [8]. Hence, in this so-called 
salt saturation limit, the bending rigidity of DNA is en- 



tirely due to the mechanical stiffness of the double-helical 
backbone. 




FIG. 1. The schematic picture of double-helical B-DNA 
with the negative charges lying on the sugar-phosphate back- 
bone in a periodic manner. 

Similar studies have shown that the twist rigidity of 
DNA is also relatively insensitive to the ionic strength 
of the solution [9]. This experimental finding is usu- 
ally explained by saying that (unlike bending) twisting a 
polyelectrolyte does not change the distance between the 
different charges on its backbone appreciably, and thus 
it is not affected by electrostatic interactions [10]. Here, 
we set out to revisit this line of argument and attempt 
to account for the above experimental observation from a 
different point of view. We consider the electrostatic self- 
interaction of the double-helical sugar-phosphate back- 
bone (see Fig. 1) and show that the periodic arrangement 
of the charge distribution effectively screens the electro- 
static interaction with the screening length given by the 
pitch of the DNA. In other words, corresponding to such 
a periodic charge distribution, there are two competing 
screening lengths: (1) the Debye screening length of the 
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bulk solution re -1 that is controlled by the ionic strength, 
and (2) the period of the charge distribution P, and it is 
the smaller of these two lengths that controls the range of 
Coulomb interaction. We find that electrostatic interac- 
tions make an appreciable contribution to the twist rigid- 
ity of DNA, although it depends only weakly on the De- 
bye screening length as long as this length is larger than 
the DNA pitch. We study the effect of various geometri- 
cal parameters such as the diameter of the double-helix, 
the distance between the two helices, and the pitch, as 
well as the Debye screening length, on the electrostatic 
contribution to twist rigidity and show that it can be 
both positive and negative depending on the values of 
these parameters. The results are summarized in Fig. 
2, where a diagram is sketched in the parameter space 
delineating all the different regimes. 



which the sugar-phosphate charged backbone of each 
DNA strand is assumed to wrap around a cylinder of 
radius a in a helical manner, as shown in Fig. 1. The 
double-helix can then be viewed as a cylinder with a sur- 
face charge density a(z, 9) corresponding to the negative 
charges, whose electrostatic self-energy can be calculated 
as 



■Eeiec = -7T / dzdz' / d6d6' a{z, 9) a(z', 9') 



x V DH (\r(z,9)~r(z\9')\) 



(1) 



where r(z, 9) parameterizes the position on the surface of 
the cylinder with z being the coordinate along the axis 
and 9 being the polar angle. The effective pair poten- 
tial between two charges in the solution is given by the 
Debye-Huckel interaction [11] 
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FIG. 2. The diagram delineating the different regimes 
in the parameter space of a helical polyelectrolyte, where 
ojq = 2n/P, and a is the radius of the helix. The line separat- 
ing the two screening regimes has slope one, whereas the slope 
of the boundary denoting the onset of instability in the Debye 
screening regime is set by the inverse of a cutoff number n c 
(see below). 

The rest of paper is organized as follows: Section II de- 
scribes the model that is used to study the electrostatic 
contribution to twist rigidity of DNA, followed by the 
presentation of the results in Sec. III. Finally, Sec. IV 
concludes the paper, while some details of the calcula- 
tions appear in three Appendices. 



II. THE MODEL 

To study the effect of electrostatic interactions on the 
twist rigidity of DNA, we consider a simple model in 



where £b — e 2 /(ek^T) is the Bjerrum length, and k 1 is 
the Debye screening length, defined via [5] 



(3) 



where Zi and Ci are the valence and the concentration of 
the salt species i, respectively, and summation holds over 
the ionic species in the solution. 

Due to the helical structure of DNA, a(z, 9) is a dou- 
bly periodic function, namely, a(z, 9) — a(z + P,9) — 
a(z, 9 + 2tt), where P is the helix pitch. Therefore, it is 
convenient to write a(z, 9) in the Fourier space as 



a(z,9) = a mn e 

m , n 

where, m and n are integer numbers 



z-\-inO 



(4) 




2 



FIG. 3. The schematic picture of the surface charge distri- 
bution of B-DNA. The geometrical parameters of DNA have 
been shown in the picture: b is the vertical distance between 
two successive charges on each strand, £ is the distance be- 
tween the two strands along the z-axis (given by the width of 
the minor groove in B-DNA), and P is the pitch of the helix. 

Making use of the periodicity of the charge distribu- 
tion, one can simplify the form of the electrostatic self- 
energy of Eq. (1) using the Fourier representation of the 
screened Coulomb interaction. After some manipulation, 
whose details can be found in Appendix A, one finds 



PE cicc = 4ir 2 £ B La 2 ^ |er m ,„| 2 /„ (\J {na) 2 + (nau) ) 

m,n ^ 

x K n (\J {no) 2 + (nau) ) 2 ) , (5) 



where (3 = 1/(/cbT), u> = 2ir/P is the spontaneous twist 
of the helix, and L is the overall length of the macro- 
molecule. 

We now focus on the specific case of DNA, whose 
charge density a(z, 9) can be written as (see Fig. 3) 



a(z,9) = - 



2nab 



(6) 



where b is the vertical distance between two charges in 
a strand, and £ is the distance between the two strands 
along the z direction, as shown in Fig. 3. Note that 



P r2ir 

dz / ad9 cr(z, 
Jo 



2P 

'T' 



(7) 



yields the number of charges in each repeat unit of DNA. 
The Fourier transform of the charge density a m ,n can 
now be calculated from Eqs. (4) and (6) as 
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2nab 



[1 + e inUoC ] 



(8) 



using which the electrostatic self-energy of the double- 
helical DNA can be calculated (from Eq. (5)) as 



/3-Bolcc 



4£bL °° ( I 

—gr~ X! + cos nLJ °0 !n ( Y ( Ka ) 2 + {nauo) 

n=0 ^ 

x K n + inauJaf^j , (9) 

where the prime indicates that the n = term should be 
counted with a prefactor of 1/2. 

To calculate the contribution to twist rigidity from the 
above Coulomb interaction, we impose an additional uni- 
form twist of £1 in the double-helix and calculate the 
change in the self-energy, i.e. fiE c \ cc (u: + ^l) — /3E e \ ec (u>o)- 
Expanding the energy change in powers of f2, we can then 
read off the electrostatic twist rigidity as 



1 d 2 (3E clcc 

elec ~ l on 2 ' 



(10) 



subject to the constraint that the relative positioning of 
the two helical strands should not alter upon deforma- 
tion, and thus we should assume that the parameter £ 
changes accordingly to (' such that (u>o + ^)C' = w oC- 



III. THE RESULTS 

Under normal physiological conditions, k s=s 1 nm _1 
and the spontaneous twist of B-DNA is uo = 1.85 nm _1 . 
Since the closed form calculation of C c i cc from Eqs. 
(9) and (10) is cumbersome, we choose to expand the 
modified Bessel functions /„ and K n to forth order in 
(na) / (nau>o) . This approximation appears to yield suf- 
ficient accuracy for the experimentally relevant range of 
parameters. 

It is convenient to use the asymptotic forms of I n (nx) 
and K n {nx) for sufficiently large n. We find that 
I n {nx)K n {nx) = ^77=5 + O(^) with S > 0, and 
observe that to a good approximation one can just use 
the relevant asymptotic forms of /„ and K n for n > 2, in 
calculating the electrostatic twist rigidity (see Appendix 
B for details). We find 

2l B a 2 

Cclcc = ^ (1 + COS^oC) 

x [/o(owo) - f2(au a )(Ka) 2 + / 4 (aw )('«a) 4 ] 
2£ B a 2 (2a 2 u; 2 -l) ^ 1 

where fo(x), f2(x), and fi(x) are functions defined in 
Appendix C. In the summation term in Eq. (11) above, 
where we have used the asymptotic forms of the Bessel 
functions, no dependence on k remains and only the 
structural parameters of DNA such as a and loq enter. 
The summation diverges as 1/n, and needs to be reg- 
ularized with a cutoff for n, which can be estimated as 
n c = 2na/t, where t is set by the thickness of each strand. 
Then, Eq. (11) can be written as 



Cclcc — 



+ 



~1p~ 



(1 + COS^oC) 



x [/o(awo) - f2(aujo){Kay 
2£ B a 2 (2a 2 a; 2 -l) 
b 2 (a 2 uj 2 + I) 5 / 2 



/ 4 (aw )H 4 



7 + In 



2 sin(w C/2) 



- (1 + cos^oC) 



(12) 



where 7 = 0.577216 is the Eulcr's constant. Note that 
the above result, as we have already mentioned, is only 
valid for k < ujq. 

Let us first evaluate the overall magnitude of the elec- 
trostatic twist rigidity, as given by Eq. (12). For B- 
DNA, we have a = 1.2 ran, ui = 1.85 ran -1 , b = 3.4 A 
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and £ = 1.13 nm [2], and the Bjerrum length is given as 
£b = 7.1 A. Using these parameters, we find Co = 174 A, 
which is relatively large. To estimate n c for B-DNA, we 
use t w 5 A, which gives n c w 15. Using these esti- 
mates and k» 1 nm -1 , Eq. (12) yields C i oc = 46 A at 
physiological salt concentration. 



0.55 
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FIG. 4. C i cc /Co as a function of aa>o- This plot corre- 
sponds to Ka = 1.2 and fo> = 2.1. 



To study the effect of various parameters, namely, the 
spontaneous twist, the diameter, and the asymmetry of 
the double-helix, as well as the salt concentration, we 
choose to work with three dimensionless parameters of 
aujQ, na, and Qujq. In Fig. 4, the behaviour of C c i C c is 
shown as a function of auo for na = 1.2 and (ivo =2.1. 
The domain for au is chosen such that the condition 
k < ujq is satisfied and Eq. (12) is valid. The plot 
shows that for sufficiently high salt concentration, the 
electrostatic torsional stiffness decreases as the sponta- 
neous twist of the double-helix increases. 
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FIG. 5. Coicc/Co as a function of acoo- This plot corre- 
sponds to kb = 0.2 and (ui = 2.1. 



For sufficiently low salt concentration, however, it ap- 
pears that the behavior is not always monotonic as shown 
in Fig. 5, where C c i C c is plotted as a function of auo 
for na — 0.2 and (u; =2.1. Interestingly, one can see 
that Ccioc can even become negative, due to the fact in 
Eq. (12) the first term becomes relatively weak for low 
salt concentrations and the second term that is dominant 
changes sign for aoJo < 1/V2- 




'0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



FIG. 6. Cdoc/G) as a function of £luq/-k. This plot corre- 
sponds to na = 1.2 and auio = 2.2. 



In Fig. 6, the dependence of C c ioc is shown on the 
asymmetry parameter Qluq for na = 1.2 and aujQ = 2.2. 
One observes that for the relatively large window of 
0.47T < (ivo < 1.6-7T, Ccioc is almost constant, and is 
thus not sensitive to the relative positioning of the two 
strands. For (iv — and 27r, a divergence sets in due 
to the fact that the charges on the two strands develop 
contacts with each other. 



0.31 
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FIG. 7. Cdoc/Co as a function of na. This plot corresponds 
to (,luo = 2.1 and acoo = 2.2. 
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Finally, the behavior of C e i ec is shown in Fig. 7 as a 
function of rea for Qloq = 2.1 and aujQ = 2.2. The de- 
pendence on the Debye screening parameter in Eq. (12) 
comes only from the first term, where the negative sign 
of the coefficient of (rea) 2 causes a dip in the plot for 
Ceiec around na w 1.3, which by chance corresponds to 
the normal physiological salt concentration. However, as 
can be seen from Fig. 7, the dependence of the electro- 
static twist rigidity on rea is extremely weak as long as 
re < luq, which is a manifestation of the fact that screen- 
ing is controlled by the periodic charge distribution and 
the effective screening length is set by the pitch P that 
is shorter than re -1 in this regime. 

It is worth saying a few words about the other limit 
where k > ujq, corresponding to high salt concentration. 
In Eq. (9), one can manifestly see that in the nth term 
in the series there is a competition between re and uloq 
to control the screening. If the salt concentration is so 
high that we have re > n c u>o, the periodic structure plays 
no role and screening is entirely controlled by the Debye 
screening in the bulk. For relatively strong Debye screen- 
ing when rea > 1, we can use the simple asymptotic forms 
of the Bessel functions and find an asymptotic expression 
for the electrostatic twist rigidity as 



Ceiec 



2£ T 



b 2 



y^(l + cosnuoC) 



n 2 [2(naLU ) 2 — (rea) 



n=l 



[M 2 



(naw ) 2 ] 5/2 ' 
(13) 



This expression can be used for the region re > u , where 
it predicts C i oc > for re < sw , and C i cc < for 
re > sivo, for a value of s w n c . 



IV. DISCUSSION 

For a polyelectrolyte with a periodic spatial charge dis- 
tribution, such as the double-helical structure of DNA, 
there are two competing mechanisms for screening the 
electrostatic self-interaction and its contribution to twist 
rigidity, namely, the Debye screening due to the free ions 
in the solution, and the structural screening caused by 
the periodic structure of the charge distribution. While 
the screening length for the former case is set by the De- 
bye length re -1 , it is set by the period P of the charge 
distribution in the latter, and the dominant mechanism 
corresponds to the one with the shorter screening length. 

It appears that the contribution of electrostatic inter- 
actions to twist rigidity can be both negative and posi- 
tive, depending on the parameters. The negative values 
for the electrostatic torsional stiffness could lead to in- 
stability in the structure of the helical polyelectrolyte, 
depending on whether the mechanical structure of the 
macromolecule can counter-balance the effect of the elec- 
trostatic instability. We have used this criterion in Fig. 
2 to summarize these different regimes in the parameter 
space. Considering that a helical polyelectrolyte seems to 



be the general structure of many stiff biopolymers (such 
as DNA and actin) it will be interesting to know which 
helical configurations can in principle lead to stable struc- 
tures, and which ones cannot. This could especially be 
important in the case of biopolymers that self-assemble 
through polymerization processes, such as actin, where 
such energetic considerations could hamper or favor the 
polymerization process. 



The twist rigidity of B-DNA is believed to be roughly 
75 nm [13], which should be thought of as the sum 
of the mechanical and the electrostatic contributions: 
C = C mec h + Coioc. Our estimate of C c i oc s=s 5 nm re- 
veals that about 7% of the twist rigidity of DNA is due 
to electrostatic interactions. This result is more or less in- 
dependent of the salt concentration, as the dependence of 
Coioc to salt concentration is very weak. For example, the 
difference between C e i ec at zero and very high salt con- 
centration is about 3.5 A. Therefore AC e i cc /C = 0.5%, 
which is very small. While this naturally explains why in 
the experiments no sensitivity on the salt concentration 
has been observed, it certainly does not mean that the 
electrostatic contribution to C is negligible. 



In the above analysis we have assumed that imposing 
a finite angle of twist between the two ends of a helical 
polyelectrolyte leads to a uniform twist. This is anal- 
ogous to the assumption of a uniform bending made by 
Odijk when calculating the electrostatic bending rigidity, 
and presumably holds true when the effective elasticity 
due to electrostatics is local. In Ref. [14], this assumption 
has been scrutinized and it has been shown that this as- 
sumption is valid provided one of these conditions hold: 
(1) the polyelectrolyte segment is long, (2) the Debye 
screening is strong, (3) the charging is weak, or (4) the 
mechanical stiffness of the polyelectrolyte is larger than 
the electrostatic contribution. We expect the same ar- 
gument to hold true for the twist rigidity as well. Since 
for the case of DNA we have shown that the mechanical 
twist rigidity is much larger than the electrostatic contri- 
bution, we can safely assume that the twist is uniform. 



In conclusion, we have studied the electrostatic con- 
tribution to twist rigidity of DNA, taking into account 
its dependence on salt concentration in the solvent. We 
have shown that there is a non-negligible electrostatic 
contribution to twist rigidity, which varies very slowly 
by changing the salt concentration in the solution. By 
changing the geometrical parameters of the helix and the 
Debye screening length, the electrostatic twist rigidity 
can change sign and become negative, implying that a 
helical structure could be a both stable as well as unsta- 
ble configuration for a helical polyelectrolyte. We finally 
note that the present analysis can be also applied to other 
biopolymers such as F-actin. 
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APPENDIX A: COULOMB ENERGY IN 
FOURIER SPACE 

Due to the periodicity of the charge distribution, it is 
convenient to calculate the electrostatic self energy in Eq. 

(1) in Fourier space. We start from the Fourier represen- 
tation of the screened Debye-Hiickel interaction in Eq. 

(2) : 

d 3 k 



f 

Jo 



[J v {x)Y dx = I v {h)K v {h), (A6) 



x 2 + h 2 

we obtain the result quoted in Eq. (5) 



APPENDIX B: ASYMPTOTIC FORMS OF THE 
BESSEL FUNCTIONS 

In this Appendix, the asymptotic forms of /„ (nx) and 
K n (nx), for large n are derived. We use the integral 
representation of these functions: 



^iec = y / / dzdz' f\ede'a(z,e)a(z',e') In{nx)= i 1 pgy f +1 dp ^ (l _ v ^-\ 

J \ ) J Jo n2 (n — I)! V 2 / J_! 



4n£ L 
k 2 + k 2 



Jk-[r(z,9)-f(z',9')}^ ^ A1 ^ 



Using the cylindrical coordinates, we can write f(z, 9) = 
zz + a(cos6* x + sin 6* y), which in conjunction with 
k = k z z + k±(cos<j) x, + sin(/> y), yields 



0E, 



clcc 



dOdO' 



+oc dk z r°° k ± dk ± 



2ir J (2tt) 2 J q 



d<t> 



47T 



k 2 + k\ + K 2 



m,n m' ,n' 

J i(2p-m+k z ^)z J [-in9+ik±a cos(0 — 9) 
g*(^ m ' — k z )z' +in' 6' — ik±a cos(0— 9') (A2) 

After integration over z, z' , and k z , we find 

k±dk± 



fiE c \ cc 



x e 
x 



„2 r 2 ^ 



dOdO'c 







m,n m' ,n' 



o (2^) 2 

47rL(5 m ,_ m / 



x e 



(n0+n'0')+iafcx[cos(>-0)-cos(4>-0')] 



(A3) 



By defining 9 = 9\ + (f> and & = 9 2 + <fi, the integration 
over the three angles in Eq. (A3) can be performed as 

J d9d9 / d(j)e i ( ne+n ' e '^ +iak± [ cos (4>-s)-cos(4>-9')] _ 

(2n) 3 5 n ,- n ,[Mk^a)] 2 , (A4) 

to yield 

f3E clcc = 4ir 2 £ B La 2 Wm.nl 2 



x / dk±k± 



J%(k±a) 



{^m) 2 + k\ + k 2 ' 



(A5) 



K n (nx) = (^) " / dp e™> (p> l)""* . 



,)! V2 



(Bl) 



Let us first consider the case of I n (nx). We denote 

r +i ! 
Q = y dp e nxp (l - p 2 ) n 1 = J dpe 9(p \ (B2) 

and expand gr(p) around po> the position of its maximum, 
to second order of (p — po), as 



9{p) 9(Po) + ^f"(Po)(p - Po) 2 
1 



= nxp + [ n - - j In (1 - p ) 

^Ap-Po) 2 , (B3) 
2 / (1-Po) 



1\ 
n 1 



where 



-l + VTT^ 2 " i / 1 :i 

= 1- - — I 1 



x 2nx \ y/l + x' 

Using this form for g(p), Q can be found as 



(B4) 



O ~ 1 Po _ nxp + (n-±)\n(l-pl) fgc^ 

Vn(l+P 2 ) 

Finally, using the Stirling's formula for nl, we find 
I n {nx) ~ — = Pa = 

x e («-i) ln ( 1 -Po)+«^Po+n(l+lnf) + ^^ (gg) 

Using a similar treatment, we find the large n asymptotic 
behavior for K n (nx), as 



G 



K n (nx) 



1 p'^-l 



2n ^n(l+p(, 2 ) 
x e("~i) ln ( p " 2 - 1 )-" 2; Po+"( 1 + ln f ) + 2i^ (B7) 



where 



, i + \/iT^ 2 i 



Po = 



2nx V vT+~ 



(B8) 



By using these relations for I n (nx) and K n (nx), we find 



I n (nx)K n (nx) 



' L +°(rir)> (B9) 



2\/l + x 2 n n 




n=2 

- - n=3 

n=5 



4.5 



FIG. 9. DifTorcnce between the asymptotic form of 
I„(nx)K n (nx) and its exact value as a function of x for dif- 
ferent values of n. The solid line corresponds to n — 2, the 
dashed line corresponds to n = 3, and the dotted line corre- 
sponds to n — 5. The difference is less than 0.3%. 

In Fig. 9, we show the behaviour of A(n, x) as a func- 
tion of x for different values of n. This plot shows that 
for n > 2, the difference goes to zero as x increases. 



where 5 > 0. 




n 

FIG. 8. Difference between the asymptotic form of 
I„(nx)K n (nx) and its exact value as a function of n for several 
values of x. The solid line corresponds to x = 0.8, the dashed 
line corresponds to x = 1.8, and the dotted line corresponds 
to x = 2.8. The difference is less than 0.25%. 

We define A(n,x) = I n (nx)K n (nx) } = , and in 

Fig. 8 show the behaviour of A(n, x) as a function of 
n for several values of x. As can be seen, the difference 
is less than 0.25% in the worst case, which implies that 
the asymptotic form of I n (nx)K n (nx) for n > 2 can be 
used as a good approximation for the range of x we are 
interested in. 



APPENDIX C: THE EXPLICIT FORMS OF THE 
AUXILIARY FUNCTIONS 

In this Appendix, we give the explicit forms of the aux- 
iliary functions used in Eqs. (11) and (12) above. The 
function fo(x) reads 

fo(x) = 4 (\ + h(x)Ki Or) +AI[{x)K[{x) 

-^[hWKiW+IiWK^x)], (CI) 

where the prime indicates differentiation. One can show 
that for x <C 1 it behaves as fo(x) = Inx + 0(x 2 lna;). 
The second function f2(x) is written as 

f 2 (x) =-3^ + ^3) lh(x)K[(x) +I[(x)K 1 (x)} 

+ - x [I[{x)K' 2 {x)-I' 2 {x)K[{x)] 

+ ^[h(x)K 1 (x)-2I[(x)K' 1 (x)] 

+ ± I [I' 2 { x )K 1 {x)-h{x)K' 2 {x)] 

+ ^[h{x)K 2 (x)-I 2 (x)K 1 {x)}, (C2) 

and for x <C 1 it behaves as f 2 (x) = l/(2.x 2 ) — lnx/2 + 
0(x 2 lna;). Finally, for f4(x) we have 

176 + 136x 2 + 19x 4 
A(a:) = ^ htyKiix) 
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-—{I^KiW+I^K^x)} 

+ ±[h(x)K' 1 (x)-I[(x)K 2 (x)] 

1 1 j_ Or 2 1 
+ —^I'MKix) + —l 3 ( x )K 3 (x) 

+ [I^K^x) - h{x)K' 2 {x)] 

88 4- 4Qr 2 

^-5 {x)Ki (x) + h {x)K[ (x)\ 



8x 5 



- — \I 2 {x)Ki{x)-h{x)K 2 {x)], 



(C3) 



which behaves as f 4 (x) = 3/(4x 4 ) - l/(8x 2 ) + 5/641nx + 
0(a; 2 ln:z!), for x -C 1. At infinity, all of these functions 
go to zero faster than 1/a; 2 , with f&{x) vanishing faster 
than f2{x), and fzix) faster than fo{x), respectively. 
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FIG. 10. The auxiliary functions fo(x), f2(x), and fi(x). 
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